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Within the effective mass approximation an adiabatic description of spheroidal 
and dumbbell quantum dot models in the regime of strong dimensional quanti- 
zation is presented using the expansion of the wave function in appropriate sets 
of single-parameter basis functions. The comparison is given and the peculiar- 
ities are considered for spectral and optical characteristics of the models with 
axially symmetric confining potentials depending on their geometric size mak- 
ing use of the total sets of exact and adiabatic quantum numbers in appropriate 
analytic approximations. 
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1 Introduction 

To analyze the geometrical, spectral and optical characteristics of quantum dots 
in the effective mass approximation and in the regime of strong dimensional 
quantization following 1J, many methods and models were used. We men- 
tion some of them, that are in the field of our interest: the exactly solvable 
models of spherical and cylindrical layer (toroid) impermeable wells [2J [3] , the 
adiabatic approximation for a lens-shaped well confined to a narrow wetting 
layer 0] , and a hemispherical impermeable well [5] , the model of strongly oblate 
or prolate ellipsoidal impermeable well [3 [8] , as well as numerical solutions of 
the boundary value problems (BVPs) with separable variables in the spheroidal 
coordinates for wells with infinite and finite wall heights [HI EH EH EH IT5] . 
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Mobius [13] nanostructures, scattering problems for toric [TS] and coupled non- 
identical microdisks [TS] . 

Similar models were used for describing the energy spectra of deformed nu- 
clei [HI HI mUni 123 121 122, atomic clusters deposited on planar surfaces [23] 
and low-energy barrier nuclear reactions [23 [23 123 123 ES] • However, thorough 
comparative analysis of spectral and optical characteristics of models with dif- 
ferent potentials, including those with non-separable variables, remains to be a 
challenging problem. 

In the present paper we analyze the spectral and optical characteristics of the 
following models: a spherical quantum dot (SQD), an oblate spheroidal quantum 
dot (OSQD), a prolate spheroidal quantum dot (PSQD), and a dumbbell QDs 
(DQD). We make use of the Kantorovich method that reduces the problem to 
a set of ordinary differential equations (ODE) [31] by means of the expansion 
of the wave function in appropriate sets of single-parameter basis functions }32j 
similar to the well-known adiabatic method [33] . 

We present briefly a calculation scheme for solving elliptical BVPs with 
axially-symmetric potentials in cylindrical coordinates (CC), spherical coordi- 
nates (SC), oblate spheroidal coordinates (OSC), and prolate spheroidal coor- 
dinates (PSC). Basing on the symbolic- numerical algorithms (SNA) developed 
for axially-symmetric potentials [531 133 ISS] > different sets of solutions are con- 
structed for the parametric BVPs related to the fast subsystem, namely, the 
eigenvalue problem solutions (the terms and the basis functions), depending 
upon the slow variable as a parameter, as well as the matrix elements, i.e., the 
integrals of the products of basis functions and their derivatives with respect to 
the parameter. These terms and matrix elements form the matrices of variable 
coefficients in the set of second-order ODE with respect to the slow variable, 
which are calculated in special cases analytically and in the general case us- 
ing the program ODPEVP [37]. The BVP for this set of ODEs is solved by 
means of the program KANTBP [38] , while in the special cases crude diagonal 
estimations can be performed using the appropriate analytic approximations. 

The efficiency of the calculation scheme and the SNA used is demonstrated 
by tracing the peculiarities of spectral and optical characteristics in the course 
of varying the ellipticity of the prolate or oblate spheroid and dumbbell in the 
models of quantum dots with different confining potentials, such as the isotropic 
and anisotropic harmonic oscillator, the spherical and spheroidal well with finite 
or infinite walls approximated by smooth short-range potentials, as well as by 
constructing the adiabatic classification of the states. 

The paper is organized as follows. In Section [2| the calculation scheme 
for solving elliptic BVPs with axially-symmetric confining potentials is briefly 
presented. Sections [3] and [4] are devoted to the analysis of the spectra and ab- 
sorption coefficient of quantum dot models with three types of axially-symmetric 
potentials, including the benchmark exactly solvable models. In Conclusion we 
summarize the results and discuss the future applications. 
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2 Problem Statement 



Within the effective mass approximation under the conditions of strong dimen- 
sional quantization, the Schrodinger equation for the slow envelope of the wave 
function 'J'(r) of a charge carrier (electron e or hole h) in the models of QDs 
has the form [5J [7] 

{H - £}*(F) = {(2^ p )- 1 P + U{f) - E}V(f = 0, (1) 

where r e R 3 is the position vector of the particle having the effective mass 

/i p = /i e (or jU p = jjihji P — — i/iVs is the momentum operator, E is the energy 
of the particle, U (r) is the axially-symmetric potential confining the particle 
motion in SQD, PSQD, or OSQD. In Model A, U (f) is chosen to be the potential 
of an isotropic or anisotropic axially-symmetric harmonic oscillator in Cartesian 
coordinates r = {x, y, z}: 

U A f) = /^ 2 (d(i 2 + f) + Csz 2 )/2. (2) 

Here £i = 1, C3 = 1 for a spherical QD or d = (f /a) 4 , £3 = (f /c) 4 for a 
spheroidal QD, inscribed into a spherical one, where a and c are the semiaxes 
of the ellipse which transforms into a sphere at a = c — Tq — \Jx q + + Zq , 
uj = "fr h/ (Pp?o) i s the angular frequency, and jf is an adjustable parameter. 
We will use the value 7^ = 7r 2 /3 that follows from equating the ground state 
energies for the spherical oscillator and the spherical QD of Model B considered 
below. If necessary, this definition can be replaced with a different one, e.g., the 
one conventional for nuclear physics [5TJ [551 [53] . 

For Model B, U(r) is the potential of a spherical or axially-symmetric well 

l> B (rO = {0, S{f) < 0; U Q , 3(f) > 0}, (3) 

bounded by the surface S(f) — with walls of finite or infinite height 1 <C Uo < 
00. In Eq. ^ 3(f depends on the parameters a, c, and < 61 < 1 

{T) - fi 2 + 5 2 (£ 2 £ 2 /4 + l-£ 2 ) 2 {) 

At ci = we get a spheroidal quantum dot model, at < c\ < 1 it becomes a 
dumbbell QD with a symmetric double well, and at c\ > 1 we get a triple-well 
model. 

For Model C, U(r) is taken to be a spherical or axially-symmetric diffuse 
potential 

U G f) =U (l-(l + cxp^r*)/*))^ , (5) 

where s is the edge diffusiveness parameter of the function smoothly approx- 
imating the vertical walls of finite height Uq. Below we restrict ourselves by 
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considering Model B with infinite walls Uq — > oo and Model C with walls of 
finite height Uq. 

Throughout the paper we make use of the reduced atomic units [TJ [7]: 
a* B = nil 2 / p, p e 2 is the reduced Bohr radius, k is the DC permittivity, Er = 
Ry* = h 2 /(2p, p a* B 2 ) is the reduced Rydberg unit of energy, and the following 

dimensionless quantities are introduced: &(r) = a* B ^ 3 ^ 2l ^(r), 2H — H/Ry*, 
2E = E/Ry*, 2U(r) = U{f)/Ry*, r= ~r/a* B , a = a/a* B , c = c/a* B ,Ci = ci/a* B , 
r a = ro/flg, oj — "ir / r Q = hui/(2Ry*). For an electron with the effective 
mass fi p = fi e — 0.067mo at n = 13.18 in GaAs: a* B = a e B — 104 A = 10.4 
nm and Ry* = E B — 5.275 meV. For a heavy hole with the effective mass 
Hh = jLt e /0.12 = 0.558mo the corresponding values are a g — a B (fi e / fXh) — 
12.48 A = 1.248 nm, and E% = E^(fj,h//j, e ) = 46.14 meV. 

Note, that for model A of approximation of OSQD/PSQD by the anisotropic 
oscillator ([2 ) the separation of variables in cylindric coordinates x = (z, p, ip) is 
possible and additional integrals exist [3j|l SOI SI] ■ Similarly, for model B the 
variables are separable in the oblate/prolate spheroidal coordinates x = (£, 77, <p) 
and the additional integrals of motion are A: [H, A] = HA — AH = 0, i.e. H p 
and A p in PSQD 
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Eq. (pi) is obtained by substituting £ — > i£, d 
in g3l3]. 

Since the Hamiltonian H in Eqs. |l])-([5| commutes with the z-parity op- 
erator of reflection in the plane z = (z — > —z or 77 — > —77), the solutions are 
divided into even (a = +1) and odd [a = —1) ones. The solution of Eq. (JlJ, pe- 
riodical with respect to the azimuthal angle ip, is sought in the form of a product 
y(xf,x s ,<p) = ^ ma (x f ,x s )e imip /V^, where m = 0,±1,±2, ... is the magnetic 
quantum number. Note, in absence of magnetic fields the Hamiltonian com- 
mutes also with the inversion operator (r — > — r) with eigenvalues a = (— 1)"V 
and solution divided into gerade (a — +1) and ungerade {& = —1) ones. Then 
the function ^ ma {xf, x s ) satisfies the following equation in the two-dimensional 
domain ft = Q Xf (x s ) U Q Xs C R 2 \{0}, il Xf (x.) = (xf a (x s ),xf a *(x s )), tl Xs = 



mm T max 
5 x i 



H^XfiXs) + H 2 {x s ) + V(x f ,x s ) - 2E) ^ ma {x fl x s ) = 0. (10) 
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Table 1: The values of conditionally fast Xf and slow x s independent variables, 
the coefficients gi s (x 8 ), gjf(%f) and the potentials Vfixf), V s (x s ), Vf s (xf,x s ), 
in Eqs.@-@ for SQD, OSQD and PSQD in cylindrical (CC), spherical (SC) 
and oblate & prolate spheroidal (OSC & PSC) coordinates with (d/2) 2 = ±(a 2 — 
c 2 ), + for OSC, - for PSC. 
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The Hamiltonian of the slow subsystem 77 2 (^ s ) is expressed as 

1 d c) 

H 2 (x s ) = H 2 (x B ) = r -r 7 7—92s(x s )-— + V s (x s ), (11) 

g ls (x s )dx s dx s 

and the Hamiltonian of the fast subsystem Hi(xf,x s ) is expressed in terms of 
the reduced Hamiltonian Hf(xf]x s ) and the weighting factor g3 S (x s ): 

Hi(x f ;x s ) = g^(x s )H f (xf,x s ), (12) 

Table [T] contains a detailed description of the conditionally fast Xf and 
slow x s independent variables, the coefficients gi a (x 3 ), g2 S (x s ), g3 S (x s ), gif(xf), 
g%f{%f), and the reduced potentials Vf(xf), V s (x s ), Vf a (xf,x s ), entering Eqs. 
([TO])-^ for SQD, OSQD, and PSQD in cylindrical (x = (z,p,<p)), spherical 
(x = (r,r) = cos6,ip)), and oblate/prolate spheroidal (x — (^,T],(p)) coordinates 
(CS, SC and OSC/PSC) 01]. Note, that in Table [TJ using Eqs. {2}, {HJ in the 
reduced atomic units, the potential V(r,Tj) for OSQD/PSQD in SC is expressed 
for Model A as 

V(r, V ) = 2r 2 U A (r, v ) = c 2 r 4 (d(l - r, 2 ) + Ca^ 2 ), 
and for Model C as 

V(r, n) = 2r 2 U c (r, „) = 2r 2 U (l - (l + exp((r 2 ( ^i^l + J) - l)/s) 

both having zero normal first derivatives dV(r, rj)/dr in the vicinity of the origin 
r = (equilibrium point), similar to [26]. We do not use the CC for Model C, 
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because the motion in this case is not restricted by two coordinates p and z. 
For Model B in Tabic [j] w = and the potentials V(r, 77) = V"(f , rj) = are 
zero, since in this case one should impose the Dirichlet boundary conditions 
^ ma (x f ,x s )\ dn = at the boundary dtt = {K 2 \S(x f ,x s ) = 0} of Q restricted 
by the surface 5(r) = 0, which is equivalent to the action of the potential (|3j) 



The solution *f l0 '(x / , x a ) ee *f m<7 (a;/, x s ) of the problem ([lOjhQ is sought 
in the form of Kantorovich expansion [31j 

^(x f ,x s )JY,^(x f ;x s ) X ^ l \E,x s ). (13) 
3=1 

The set of appropriate trial functions is chosen as the set of eigenfunctions 
$y i<T (x/; x s ) of the Hamiltonian Hf{xf\x s ) from (12 1, i.e., the solutions of the 
parametric BVP 

{H f (x f ; x.) - k(x s )} $r = 0, (14) 

in the interval Xf £ X/ (a; s ), depending on the conditionally slow variable x s € 
Q Xa as on a parameter. These solutions obey the boundary conditions 

lim /iVf^g^ (15) 
at the boundary points {x™ n (:E s ), x™ ax (a; s )} = (30^ (# s ), of the interval (x s ). 



In Eq. jl5|, 7V} m<T) (x s ) ee ATj m<T) , D*"^^) ee £>J m<T) , unless specially declared, 



are determined by the relations N^" 1 ^ = 1, = at m = 0, a = +1 (or 

at er = 0, i.e., without parity separation), 7vj mCT ' > = 0, D^"" 7 ' = 1 at m = 0, 
<t = — 1 or at m 7^ 0. The eigenfunctions satisfy the orthonormality condition 
with the weighting function gifixf) in the same interval x/ € f2 X/ (x s ): 

•/z™ ln (x s ) 

Here \i(x s ) < ... < A Jmax (a; s ) < ... is the desired set of real eigenvalues. The 
corresponding set of potential curves 2E 1 (x s ) < ... < 2£ , Jmax (x s ) < ... of Eqs. 



Q is determined by 2E ] (x s ) = g^ 1 (x s )\j(x s ). Note that for OSC and PSC, 
the desired set of real eigenvalues Xj(x s ) depends on the combined parameter, 
x s — > p 2 = (d/2) 2 2E, i.e., the product of spectral 2E and geometrical (d/2) 2 
parameters of the problem (10 1. The solutions of the problem (14|-(fl6) for 
Models A and B are calculated in the analytical form [36] , while for Model C 
this is done using the program ODPEVP j37j. Substituting the expansion (13 1 
into Eq. ([IJ , we get a set of ODEs for the slow subsystem with respect to the 

unknown vector functions x^ mat ^ (x s , E) ee x^\x s ) = [x\ ( x s), —iXj ( x s)) T '- 

{ 1 1 —^g2s{x s )^ + Vs(x s ) + Vu(x s )-2E\ x ( i ) (xs) =-Y J Vij{x a ) x f(x4l7) 

\ gis(Xs) dx s dx a j j 
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Here Vu(x s ) = 2Ei(x s ) + Ha(x s ), Vij(x s ) are defined by formula 

gi s {x 3 ) gis(xs) dx s gis(xs) dx s 



H lj {x a ) = H jl {x s )= g xf {x f ) ,( J' 81 3 \1' s, dx f , (18) 

J x l j? la (x s ) 



dx s dx a 



Qij(xs) = -Qji(x a ) = — I _" g 1 f(xf)$i(xf,x 3 ) J — — - dxj 

(x B ) OXs 



r x f 1 

J x J ^ LITl (x 

and calculated analytically for Model B and by means of the program ODPEVP 
[37] for Model C, while the solutions of the BVPs for Eq s. 0\ with the bound- 
ary and orthonormalization conditions of the type ( 15 ) , ( 16 ) with Xf — > x s were 
calculated by means of the program KANTBP |38). Note, that for Model A 
in SC or CC and Model B in OSC or PSC, the variables Xf and x s are sepa- 
rated so that the matrix elements Vij(x s )0 are put into the r.h.s. of Eq. (17), 
and V s (x s ) are substituted from Table [T] For the interesting lower part of the 
spectrum of Models A and B 2E : 2E 1 < 2E 2 < . . . < 2E U or of Model C 
2E : 2Ei < 2E 2 < . . . < 2E t < 2Uq, the number j max of the equations solved 
should be at least not less than the number of the energy levels of the problem 
( 17| at a = c = ro • To ensure the prescribed accuracy of calculation of the 
lower part of the spectrum discussed below with eight significant digits we used 
imax = 16 basis functions in the expansion (8) and the discrete approximation of 
the desired solution by Lagrange finite elements of the fourth order with respect 
to the grid pitch ^ s(a , s) = [x s - min ,x s . k = x s . k _i + h s ,x s . max ]. The details of 
the corresponding computational scheme are given in |36j . 



3 Spectral Characteristics of Spheroidal 
and Dumbbell QDs 

3.1 Model A of OSQD & PSQD 

In the exactly solvable model A the variables are separable in spherical co- 
ordinates, and under the variation of the aspect ratio parameters Cca = c/a 
and ( ac = — a/c for the oblate and prolate spheroids, determining the 
transverse w p = VCi^ an d longitudinal uj z — y/C,^ frequencies of the circular 
and linear harmonic oscillators. The spectrum is given by the sum of energies 
2E npTn = 2uj p (2n p + \m\ + 1), n p — 0, 1, . . . , m — 0, ±1, . . . (with the eigenval- 
ues being degenerate with respect to X p = 2n p + \m\ that number in ascending 
order the energy values of the states [45j ES] that is conventionally used in prac- 
tice, see, for example, [HI [23) an d 2i?„ z = 2uj z (n z + 1/2), n z — 0,1,..., at 
u = (j ro = 7r 2 /(3r2), ^ 1 = r 2/(a 2 ), and ^C 3 = r o/( c2 )- At a = c = r o the 
independent variables are separable in the boundary problem for Eq. ([T]) in the 
spherical coordinates too, i.e., we have the energy spectrum of a spherical oscil- 
lator 2E° s r c lm = 2uj ra (2n r + 1 + 3/2), n r = 0, 1, . . I = 0, 1, . . ., m = 0, ±1, . . . , ±1 
with the eigenvalues being degenerate with respect not only to to, but also to 
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1 1.5 2 2.5 1 1.5 2 2.5 
c a 

a) b) 

Figure 1: Energies 2E — E/E R of the even u = +1 lower states of Model A 
OSQD at a = 2.5 (a) and PSQD at c = 2.5 (b) versus c or a. The exact intersec- 
tions of the energy levels take place at rational ratios R = oj p /oj z = (c/a) 2 6 Q 
(a) and R — uj z /uj p = (a/c) 2 e Q (b) of the frequencies of transverse and 
longitudinal oscillators with R = 1,4/5,3/4,2/3,3/5,1/2,2/5,1/3,1/4,1/5,... 



A,. = 2n r + I that number in ascending order the energy values of states, sepa- 
rated in parity a = (-1) A = (-1)' = (-l)' m cr, cr = (-1)'"™ = ±1. The energy 
spectrum of the spherical oscillator 2E° sc lm coincides at a = c with 

2E(a,c) = 2(E„ eo +E npom ) and 2E(c, a) = 2(E npp , m + E nzp ), (19) 

which, respectively, defines the one-to-one correspondence between the sets of 
the quantum numbers n zo = I — \m\, n po = n r , m — m for OSQD and SQD 
and = I — \m\ for PSQD and SQD, that characterize the 

fast and slow subsystems at continuous variation of the parameters £ ca = c/a 
and C ac = a/c. At decreasing the parameter £ ca or (^ ac the degeneracy of the 
spectrum with respect to the quantum numbers n, I, m is removed. 

Fig. [IJillust rates the lower part of the equidistant energy spectrum E/Er = 
2E(a, c) and E/Er = 2E(c, a) for even states a = +1 of the model of OSQD and 
PSQD with parabolic confining potentials |2|, at m = 0, i.e., of an oblate and 
prolate spheroid, depending on the minor c or a and the major a or c semiaxes, 
respectively. At fixed values of the parity a and the magnetic quantum number 
m when the ratio of the frequencies w p and w z of the longitudinal and transverse 
oscillators is a rational number, uj p /uj z £ Q, as illustrated, e.g., in Fig. [TJ the 
exact crossings of the same-parity terms occur, after which above each energy 
level of OSQD (or PSQD), labelled with the quantum number n zo (or n pp ) 
of the fast subsystem, an equidistant spectrum appears with the energy levels 
labelled with the quantum number n po (or n zp ) of the slow subsystem. Note, 
that when the parameters tend to zero, the longitudinal energy of OSQD and 
the transverse energy of PSQD tend to infinity. However, since the variables 
are separable and the energy can be presented as a sum, the finite energies for a 



disc E nporn or a wire E nzP result from the subtraction of the longitudinal E nzo 
or transverse E nppTn energy, respectively. 

3.2 Models B and C for Oblate Spheroidal QD. 

At a fixed coordinate x s of the slow subsystem, the motion of the particle in 
the fast degree of freedom Xf is localized within the potential well having the 

effective width 

L(x s ) = 2c^l-x 2 Ja\ (20) 



where L = L/a* B . The parametric BVP for Eq. (12 1 at fixed values of the 
coordinate x s , x s £ (0, a), is solved in the interval Xf £ (—L(x s ) /2,L{x s ) /2) 
for Model C using the program ODPEVP, and for Model B the eigenvalues 
En a (%s) /En = 2Ei (x s ), n = i = 1, 2, and the corresponding parametric 
eigenfunctions $f (xf;x s ), are expressed in the analytical form: 



2E ^LW7y ^^^^V^'H^Im^" 1 ))' (21) 

where the even solutions a = +1 are labelled with odd n a = n zo + l = 2i — l, and 
the odd ones a = —1 with even n — n zo + 1 = 2i, i = 1, 2, 3, ... . The effective 



potentials (18 1 in Eq. (17) for the slow subsystem are expressed analytically in 



terms of the integrals over the fast variable Xf of the basis functions (21) and 
their derivatives with respect to the parameter x s including the states of both 
parities a = ±1: 

a 2 7r 2 u 2 3 + 7T 2 7l 2 a; 2 

2Ei(x s ) — — — - — , Ha(x s ) — — — ^ »Tj, (22) 
4c*(a /i — xj) 12 \a l — x%) £ 



2n n' (n 2 . + n' 2 )(l + (-l)"°+"°) x\ 

(n 2 — n' D 2 ) 2 (a 2 — x 2 ) 2 ' 



VvW'i (n 2 -n' 2 ) 2 a 2 -x 2 ' n "^ n °- 



For Model B at c = a = vq the OSQD turns into SQD with known analyti- 
cally expressed energy levels E t = E^ m and the corresponding eigenfunctions 

2E s P _ <+i.i+i/2 SP , Q v V^J t+ i /2 (V2^ m r) 

(23) 

where asn r +i,l+:i./2 are zeros of the Bessel function of semi-integer index I + 1/2, 
numbered in ascending order < an < ai2 < ... < < ... by the integer 
i, v = 1, 2, 3, .... Otherwise one can use equivalent pairs iv -H> {n r , 1} with n r = 
0,1,2,... numbering the zeros of the Bessel function and / = 0,1,2,..., being 
the orbital quantum number that determines the parity of states a = (—1)' = 
(-l) m cr, a = (-l) l -' m = ±1. At fixed I, the energy levels E nlm /E R = 2E t 
degenerate with respect to the magnetic quantum number m, are labelled with 
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Figure 2: The energies 2E = E/Er of even a = +1 lower states for OSQD 
versus the minor c, £ ca = c/a e (1/5, 1) being the spheroid aspect ratio: a) well 
with impermeable walls, b) diffusion potential with 2Uq — 36, s — 0.1, the major 
semiaxis a = 2.5 and m — 0. Tine lines are minimal values 2Ef im = 2Ei(x s = 0) 
of potential curves. 



the quantum number n — n r + l = i = 1, 2, 3, ... , in contrast to the spectrum of a 
spherical oscillator, degenerate with respect to the quantum number A = 2n r +l. 
Figs. [i]|3]show the lower part of the non-equidistant spectrum E(( ca ) /En = 2E t 
and the eigenfunctions <J/™ cr from Eq. ( 13 ) for even states OSQD Models B and 
C at m = 0. There is a one-to-one correspondence rule n Q — n zo + 1 = 2n— (1 + 
cr)/2, n = 1, 2, 3, between the sets of spherical quantum numbers (n, I, m, a) of 
SQD with radius tq = a = c and spheroidal ones {n^ = n r , n v — I — \m\, m, cr} 
of OSQD with the major a and the minor c semiaxes, and the adiabatic set 
of cylindrical quantum numbers [n zo , n po , m, a] at continuous variation of the 
parameter C ca = c/a. The presence of crossing points of the energy levels of 
similar parity under the symmetry change from spherical Cca = 1 to axial, i.e., 
under the variation of the parameter < C ca < 1, in the BVP with two variables 
at fixed m for Model B is caused by the possibility of variable separation for Eq. 
([8| in the OSC [33] . i.e., the r.h.s. of Eq. ( 17 ) equals zero, and by the existence 
of the integral of motion Q . The transformation of the eigenfunctions occurring 
in the course of a transition through the crossing points (marked by circles) in 
Fig. [2] is shown in Fig. [3] for model B (marked by arrows) and similar for model 
C. From the comparison of these Figures one can see that if the eigenfunctions 
are ordered in accordance with the increasing eigenvalues of the BVPs, then for 
both Models B and C, the number of nodes [47] is invariant under the variation 
of the parameter c from c = a = 2.5 to c = 0.5 of the potentials (|3| and (JSJ). For 
Model B, such a behavior follows from the fact of separation of variables of the 
BVP with the potential ^ in the OSC, while for Model C further investigation 
is needed, because the coordinate system, in which the variables of the BVP 
with the potential ^ are separable, is unknown. So, at small values of the 
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Figure 3: Contour lines of the first five even-parity wave functions a = +1 in 
the xz plane of Model B of OSQD for the major semiaxis a = 2.5 and different 
values of the minor semiaxis c (£ CQ = c/a G (1/5, 1)) 



deformation parameter (£ ca for OSQD or £ ac for PSQD) there are nodes only 
along the corresponding major semiaxis. For Model C at each value of the 
parameter a there is a finite number of discrete energy levels limited by the 
value 2Uq of the well walls height. As shown in Fig. [2J5, the number of levels 
of OSQD, equal to that of SQD at a = c = ro, is reduced with the decrease of 
the parameter c (or ( C a), in contrast to Models A and B that have countable 
spectra, and avoided crossings appear just below the threshold. 

3.3 Models B and C for Prolate Spheroidal QD. 

In contrast to OSQD, for PSQD at fixed coordinate x s of the slow subsystem 
the motion of the particle in the fast degree of freedom Xf is confined to a 2D 
potential well with the effective variable radius 

Po (x s ;a,c) = a^/l - x 2 Jc 2 , (24) 
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where po (x s ) = po (x s ) /a* B . The parametric BVP for Eq. (12) at fixed values 
of the coordinate x s from the interval x s G (— c, c) is solved in the interval 
Xf € (0, po(x s )) for Model C using the program ODPEVP, while for Model 
B the eigenvalues E 1lpp+ i (x s ) /En = 2Ei (x s ), n pp + 1 = i = 1,2, and the 
corresponding parametric basis functions (j)™ ^ (xf;x s ) = $f ( x fi x s) without 
parity separation are expressed in the analytical form: 



y/2 J \rn\(^j2E npp + lt \ m \ (x s )xf) 



2 

2Et (x s ) = Q " p 9 p + 1 '' m| , $™ (x s ) = -^j^^^-r 7 " ' — n— . ( 25 ) 

where a npp +i,|m| = ^|mf +1 are positive zeros of the Bessel function of the first 
kind J| m | (xf), labeled in the ascending order with the quantum number n pp +l = 
i = l,2, .... 



The effective potentials ( 18 ) in Eq.( 17 1 for the slow subsystem are calculated 



numerically in quadratures via the integrals over the fast variable x t of the basis 



functions) 25 ) and their derivatives with respect to the parameter x s , and at 



m = may be presented in the analytical form: 



Pa(Xs) Jo Ji{ J o) Ji( J o) 



Figures|4j [5jillustrate the lower part of the non-equidistant spectrum E(£ ac ) / Er 
= 2E t and the eigenfunctions ^™ cr from Eq. (Tl3j) of even states of PSQD Models 
B and C. 

A one-to-one correspondence rule n pp + 1 = n p = i = n = n r + 1, i = 1,2,... 
and n zp = I — \m\ holds between the quantum numbers (n, I, m, a) of SQD 



n r , n 



with the radius r = a = c, the spheroidal quantum numbers {n^ 
I — \m\,m, cr} of PSQD with the major c and the minor a semiaxes, and the 
adiabatic set of quantum numbers [n p = n pp + 1, n zp , m, u] under the continuous 
variation of the parameter £ oc = a/a The presence of crossing points of similar- 
parity energy levels in Fig. [4] under the change of symmetry from spherical 
Cac = 1 to axial, i.e., under the variation of the parameter < C, ac < 1, in the 
BVP with two variables at fixed m for Model B is caused by the possibility 
of variable separation for Eq. M in the PSC [g], i.e., r.h.s. of Eq. p7| 
equals zero, and by the existence of the additional integral of motion Q. For 
Model C, at each value of the parameter c there is also only a finite number 
of discrete energy levels limited by the value 2Uq of the well walls height. As 
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Figure 4: The energies 2E — E/Er of even a = +1 lowest states for PSQD 
depending on the minor semiaxis a (£ ac = a/c 6 (1/5, 1) is the spheroid aspect 
ratio): a) well with impermeable walls, b) diffusion potential, 2Uq = 36, s = 0.1, 
for the major semiaxis c = 2.5 and m = 0. Tine lines are minimal values 
2E rnin = 2Ei{x s = 0) of potential curves 

shown in Fig. [4Jd, the number of energy levels of PSQD, equal to that of SQD at 
a = c = ro, which is determined by the product of mass p, e of the particle, the 
well depth Uq, and the square of the radius Tq, is reduced with the decrease of 
the parameter a (or £ QC ) because of the promotion of the potential curve (lower 
bound) into the continuous spectrum, in contrast to Models A and B having 
countable spectra. Note, that the spectrum of Model C for PSQD or OSQD 
should approach that of Model B with the growth of the walls height Uq of the 
spheroidal well. However, at critical values of the ellipsoid aspect ratio it is 
shown that in the effective mass approximation, both the terms (lower bound) 
and the discrete energy eigenvalues in models of the B type are shifter towards 
the continuum. Therefore, when approaching the critical aspect ratio values, it 
is necessary to use such models, as the lens-shaped self-assembled QDs with a 
quantum well confined to a narrow wetting layer [3] , or, if the minor semiaxis 
becomes comparable with the lattice constant, to proceed to models beyond the 
effective mass approximation (see, e.g. [48]). 

3.4 Models B for dumbbell QD 

For DQD at the fixed coordinate x s of the slow subsystem the motion of the 
particle in the fast degree of freedom Xf is confined to a 2D potential double 
well at < ci < 1 with the effective variable radius 

q ~\~ 1 ■ 

po (x s ) = po (x s ;a,c,ci) = -^c 2 - x 2 s 2 S ) \. (27) 

c c^c z /4 + 1 — c\ 

Fig. [6] illustrates the transformation of the prolate spheroidal shape of QD with 
c = 2.5 and a = 0.5, considered in the previous Section, into a "dumbbell" -type 
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a-0.5 a-1.0 a-1.5 a=2.0 a-2.5 



Figure 5: Contour lines of the first five even-parity wave functions a = +1 in 
the xz plane of Model B of PSQD for the major semiaxis c = 2.5 and different 
values of the minor semiaxis a (£ ac = a/c G (1/5, 1)) 

shape and the corresponding evolution of the lower part of the countable spec- 
trum E(( ac — l/5,Ci)/En — 2E t of Model B versus the deformation parameter 
Ci at a few fixed values C\ = 0, 0.25, 1 from the interval < C\ < 1. At ci — 
the discrete spectrum states are characterized by a set of exact spheroidal or 
adiabatic cylindrical quantum numbers, {n^,n^, m, <r} or [n pp , n zp , to, er]. Typi- 
cally, one can see exact crossing of energy levels having different parity (a = ±1) 
with the growth of the deformation parameter c±, which leads, first, to the 
quasidegeneracy of these energy levels and then to their exact degeneracy at 
the critical value Cx = 1. On the other hand, for small values of the defor- 
mation parameter c\ one observes, first, exact crossings (labelled with circles 
like in Figjiji above) of similar-parity energy levels, replaced with the avoided 
crossings (labelled with squares) for greater values of the deformation param- 
eter approaching the critical value c\ = 1. A similar picture was observed in 
the example of a 2D-Sinai billiard [35], a 2D-quantum billiard with the shape 
x 2 + y 2 + ex 3 = 1 and the deformation parameter e > 0, provided the so-called 
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Figure 6: a. The profile in plane z,p of closed surface generated by rotating 
of continuous curve po(z; a, c, cy) from Q about z-axis for c = 2.5, a — 0.5 vs 
Cy = 0,0.25,0.5,0.75,1. b. The energy levels of the even and odd states of 
DSQD for model B generated by Q for c = 2.5, a — 0.5 vs cy classified at 
cy = by adiabatic quantum numbers n pp , n zp and m = of PSQD. 



whispering gallery modes and considered in [501 151) , as well as in the unidirec- 
tional far-field emission of coupled nonidentical microdisks [IB] , 

In Fig. [7J we show the evolution of the first five eigenfunctions with the 
increasing deformation parameter values cy = cy = 0, 0.11, 0.99. 

The transformation of eigenfunctions when passing the avoided crossing 
points (labelled with squares) in Fig. [6] b, is shown in Fig. [7] for model B 
of DQD (labelled with arrows). Comparing these Figures, one can see that if 
the eigenfunctions are ordered in accordance with the increasing eigenvalues of 
the BVPs, then the number of nodes is not invariant under the variation of the 
parameter cy from cy = to cy = 1 in the potentials (27). In particular, in Fig. [jj 
one can see that the eigenfunction of the state [n pp = 0, n zp = 6, m = 0, a = +1] 
at Cy = 0.99 has the same number of nodes as the eigenfunction of the state 
l,n zp = 0,m = 0, a = +1] at cy = 0. Above we could already observe 



ri 



pp — x i ' l zp 

this in Fig{5]at a = 1 (up-going arrow) after several exact and avoided crossings 
of the corresponding energy levels in Fig. [frj}). At the same time, the eigenfunc- 
tion of the state [n pp = 7 n zp = 8, m = 0, a — +1] at cy = 0.99 after avoided 
crossing of the corresponding energy levels in Fig. [6)3) has the same number of 
nodes as the eigenfunction of the state [n pp = 6,n zp = 0, m = 0, <r = +1] at 
cy = 0. 



4 Absorption Coefficient for an Ensemble of 
QDs 

One can use the mentioned differences in the energy spectra to verify the con- 
sidered models of QDs by calculating the absorption coefficient K(oj ph , a, c, ) of 
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Figure 7: Contour lines of the first five eigenfunctions of model B of DSQD at 
a = 2.5, c = 0.5 and several values of c\. Light and dark inflections are positive 
and negative values of eigenfunctions and lines are eigenfunction nodes. The 
adiabatic cylindrical quantum numbers [n pp , n zp ] are given at m — 0. Crossing 
arrows mean the transformation of nodes of pair of eigenfunctions after passing 
value of parameter in which avoided crossing of corresponding pair of eigenvalues 
was taking place in Fig. [6jb. 

an ensemble of identical semiconductor QDs 52J : 

£(u;^a,c)=^it x (£^ (28) 

ty = \ J K(ha,c,)^((ha,c,))df\ 2 , W vu . = E g + E e v {a,c) + E^,(a,c), 

where A is proportional to the square of the matrix element in the Bloch decom- 
position, ^> e v {u) and are the eigenfunctions of an electron (e) and a heavy 
hole (h), and E^, are the energy eigenvalues for an electron (e) and a heavy 
hole (h), depending on the semiaxis size c, 5 for OSQD (or a,c for PSQD) and 
the adiabatic set of quantum numbers v — [n ZOl n pol m] and v' = [n' zo ,n po i , m'] 
[y = [npp, n zp , m] and v' = [n', n' zp , m']), where m! — —m, E g is the band gap 
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width in the bulk semiconductor, £j ph is the incident light frequency, W vu > is 
the inter-band transition energy for which K(ui ph ) has the maximal value. We 
rewrite the expression ( 28 1 using dimensionless quantities in reduced atomic 
units 

K(u; ph ,a,c) = I^^^.^/.^h)], 



fuy(u) = Ax - (2E g )- 1 (2E e u (a,c) + 2£&(a,c)( Wl /p e )), 

where the parameter will u be defined below, Ai = (huj ph ~ E g ) / E g is the energy 
of the optical interband transitions scaled to E g , 2E g = E g /E R = 1.43/(5.27 • 
1CP 3 ) is the dimensionless band gap width. For both electron and hole carri- 
ers the dimensionless energies 2E% — El/E e R and 2E^(fi h /fj. e ) = E^/E e R are 
expressed in the same reduced atomic units E R . 

Now consider an ensemble of OSQDs (or PSQDs) with different values of the 
minor semiaxis c = u c (or a = u p d) determined by the random parameter u = 
u a (or u — Up). The corresponding minor semiaxis mean value is c at fixed major 
semiaxis a (or a at fixed major semiaxis c) and the appropriate distribution 
function is P{u Q ) (or P(u p )). Conventionally, they use the normalized Lifshits- 
Slezov P(u) = P LS (u) [S3] or Gaussian P(u) = P G (u) distribution functions 
(/ P{u)du = J uP{u)du = 1): 



[ 3 4 en 2 exp(-l/(l-2u/3)) 
P LS ( U ) ■= < 25/3(u+3) 7 /3(3/2-„)ii/3 • 

0, 



ue (0,3/2); 
otherwise 



P G (u) := l/V2%/aexp(-(u - l) 2 /(2a 2 )), 

where u = J uP G (u)du = 1 is the mean value of u and a 2 — (J(u — u) 2 P G (u)du) 
is the variance. The absorption coefficient of an ensemble of semiconductor QDs 
with different dimensions of minor semiaxes is then expressed as 

K°(u! ph ,a,c) = I K{uj ph 7 a,c,u )P{u )du , 



K p {uj ph 1 d,d) = / K{uj ph ,a 1 ~c,Up)P(up)du p 



Taking the known properties of the (5-function into account, we arrive at the 
analytical expression for the the absorption coefficient K(uj ph ,d,c) of a system 
of semiconductor QDs with a distribution of minor semiaxes: 



k(u ph ) 



E 



I'll) 



= 1, 



df v y (u) 




du 


U=U S 



P(U S ) 



(29) 

where Kg = A 1 E g is the normalization factor, u s are the roots of the equation 
Uy{u s ) = 0. 

In particular, for Model B of OSQD or PSQD we have the interband overlap 



for OSQD, 



17 



Iu,v> — {■h + \m\{ a n pp + l,\m\) I Jl-\m\{ a n pp + l,\m\)) ,n' ip $n pp ,n' pp S m -m' f° r PSQD, 

and the selection rules n zo — n' zo , n po = n' po , and m = —m! or n pp = n' , 
n zp = n' zp and m — —m', respectively. Note that the contributions of non- 
diagonal matrix elements to the energy values are about 1% for OSQD and 
PSQD of Model B; then in the Born-Oppenheimer approximation of the order 
bmax for the absorption coefficient we get 

Uy{u) = X 1 - g£0V- 2 . (30) 

3=0 

Here the coefficients E^' are defined by 

or m = (2E g r 1 E^ul-l(a)(l + njn h ). (31) 
w P; n (c) = nn /(ac), uj z . npp (a) = a npp+lt \ m \/(ac). (32) 

E^=a 2 /4, E £> =(2n po +H + l), 

E ?o = (6n po |m|+2+6n po +6n2 o + |m| 2 + 3|m|)a- 2 , (33) 

EfJ = 3 ( 6n po + 3 1 m \ + 2 + 1 m | 2 + 6n 2 p + 6n po | m | + 4n 3 po 

+ Q\m\n 2 po + 2\m\ 2 n po )a A /2, (34) 

E?J = c 2 , £g> = (2n zp +l), fig> = +3(2n zp + 2n 2 p + l) C - 2 /4, 

^ =3(3n 2 p + 7n, p + 2n^+3)c- 4 /16. 

The coefficients of the order 6 max > 4 are calculated by the perturbation theory 
algorithms |3U [35] using exact solutions of 2D and ID oscillators with adiabatic 



frequencies w p; „ o (c) and uj z;n (a) from (32) that distinguish from conventional 
ones, for example, u> p and u z using in section [3] 1 or in |18 [ I24 j . The accuracy 
of such approximations up to i> max = 5 is about 4-6 decimal digits in compar- 
ison with the numerical results of the crude diagonal adiabatic approximation 
(CDAA) of Eqs.(17l without Hu{x s ) for the states from Fig. [2^, at c = 0.5 and 
Fig. [4^, at a — 0.5. In the case a = c = 1 the accuracy is only about two decimal 
digits in comparison with the CDAA of the exact spectrum Eq. ( [23] ) of model 
B of SQDs [52]. 

Note that in model B 2E io and 2E ip monotonically depend upon the param- 
eter u and, therefore, the algebraic equation f v v i (u) = has the only solution in 
the considered domain of definition. Using the notations X[ — Ai for b max — 1 

(2) (2) 

and X[ = Ai — E) lo , or X\ = X\ — E ip for b rnax > 2, we rewrite this equation in 
the Born-Oppenheimer approximations up to the third order b max < 3 

/„,„(«) - Ai - E^u- 2 - E^u- 1 - E^u = 0, 
which has the required roots u\ — u^ bmax ' ) : 

u[^ 2) = (2X' 1 )- 1 (E^ + ((E^) 2 + 4A' 1 ^°)) 1 / 2 )), 
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Figure 8: Absorption coefficient K/Kq from (29) consists of sum of the first 
partial contributions vs the energy A = Ai of the optic interband transitions for 
the Lifshits-Slezov distribution in first, second and third (from left to right) Born 
Oppenheimer approximations: (top panels) for assemble of OSQDs c = 0.5, 
a = 2.5 (summation by n — 1,2,3, n po = 0,1,2,3,4,5, m — 0), (bottom 
panels) for assemble of PSQDs a — 0.5, c = 2.5 (summation by n p = 1,2,3, 



0,1,2,3,4,5, m = 0). 



For the Lifshits-Slezov distribution Fig. [HJdisplays the total absorption coef- 
ficients K{u} ph ) / Kq and the partial absorption coefficients K v ^(uj ph )/Ko, that 
form the corresponding partial sum ( |29[ ) over a fixed set of quantum numbers 
v at m = —m! = 0. One can see that the summation over the quantum num- 
bers n zo (or n pp ) numerating the nodes of the wave function with respect to 
the fast variable gives the corresponding main maxima of the total absorption 
coefficients for the ensemble of QDs with distributed dimensions of minor semi- 
axis, while the summation over the quantum number n po (or n zp ) that label 
the nodes of the wave function with respect to the slow variable leads to the 
increase of amplitudes of these maxima and to appearing secondary maxima in 
the case of sparer energy levels of Model B OSQDs (or PSQDs) 

In the regime of strong dimensional quantization the frequencies of the in- 



terband transitions between the levels 



= l,n po = 0,m = for OSQD or 



n p = l,n zp = 0,rn = for PSQD in the BOl, at the fixed values a = 2.5a e 
and c = 0.5a e for OSQD or a = 0.5a e and c = 2.5a e for PSQD, are equal to 
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n o =3,n zo =0,l,. 



n o =2,n zo =0,l,... 



n o =l,n zo =0,l,... 



n„=2,n D „=0,l,... 



n' D =l,n ' zo =0,l,...n' o =2,n ' zo =0,l,... n' =3, n' zo =0,l, 



Figure 9: Schematic plots of interband transition corresponded to Fig. [HJ 



u>$ = 2.17 x 10 13 s" 1 or = 3.32 x 10 13 s" 1 (Q^ = ft-^ioo.ioo with the 
accuracy to 3% and 0.5%, respectively), corresponding to the infrared spectral 
region [SJ[7]. With decreasing semiaxis the threshold energy increases, because 
the "effective" band gap width increases, which is a consequence of the enhance- 
ment of dimensional quantization. Therefore, the above frequency is greater for 
PSQD than for OSQD, because the SQ implemented in two direction of the 
plane (x,y) is effectively greater than that in the direction of the z axis solely 
at similar values of semiaxes. Higher-accuracy calculations reveal an essential 
difference in the frequency behavior of the absorption coefficient for interband 
transitions (see Fig. |9| in systems of semiconductor OSQDs or PSQDs having 
a distribution of minor semiaxes, which can be used to verify the above models. 



5 Conclusions 

The presented examples of the analysis of energy spectra of SQD, OSQD, PSQD, 
and DQD models with three types of axially symmetric potentials demonstrate 
the efficiency of the developed computational scheme and SNA. Only Model 
A (anisotropic harmonic oscillator potential) is shown to have an equidistant 
spectrum, while Models B and C (wells with infinite and finite walls height) 
possess non-equidistant spectra. In Model C, there is a finite number of energy 
levels. This number becomes smaller as the parameter a or c (Coc or £ CQ ) is 
reduced because the potential curve (lower bound) moves into the continuum. 
Models A and B have countable discrete spectra. This difference in spectra 
allows verification of SQD, OSQD, and PSQD models using the experimental 
data [2], e.g., photoabsorption, from which not only the energy level spacing, 
but also the mean geometric dimensions of QD may be derived [BJ [TQl [IT] . The 
considered examples of calculating the absorption coefficient for ensembles of 
OSQDs or PSQD's with random minor semi-axes in model B have proved the 
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possibility of a similar verification. It is shown that there are critical values 
of the ellipsoid aspect ratio, at which in the approximation of effective mass 
the discrete spectrum of the models with finite-wall potentials turns into a 
continuous one. Hence, using the experimental data, it is possible to verify 
different QD models like the lens-shaped self-assembled QDs with a quantum 
well confined to a narrow wetting layer 4 J , or to determine the validity domain 
of the effective mass approximation, if a minor semiaxis becomes comparable 
with the lattice constant and to proceed opportunely to more adequate models 
such as [48] . 

Further development of the method, symbolic-numerical algorithms, and the 
software package is planned for solving the quasi-2D and quasi- ID BVPs with 
both discrete and continuous spectrum, which are necessary for calculating the 
optical transition rates, channeling and transport characteristics in the models 
like quantum wells or quantum wires and low-energy barrier nuclear reactions. 
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